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FOREWORD 


This report is one of several to be published from research 
conducted under NASA Contract No, NAS8-28590, entitled, "Feasibility 
of Remote Evaporation and Precipitation Estimates*" A number of 
approaches have been and continue to be followed in the conduct of 
the research. The results presented in this report represent only 
a portion of the total research effort. Other reports will be pub- 
lished as the research progresses. 

The author extends sincere gratitude to Mr. Robert Turner 
of the Aerospace Environment Division, NASA-Marshall Space Flight 
Center, under whose guidance this research was done. His encourage- 
ment and interest were essential to the undertaking and completion 
of this report. Sincere appreciation is expressed to the personnel 
of the Aerospace Environment Division for their review of the report 
and their helpful suggestions. 
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INTRODUCTION 


Broadly, changes in the constituents of air either by addition to 
and/or subtraction from them which affect the physical and/or chemical 
properties of the air, and subsequently cause detectable deterioration 
of air quality may be defined as air pollution. Pollutants include any 
natural and artificial (man-made) airborne contaminants capable of 
altering air properties. Air pollutants can occur in form of gases, 
liquid droplets and solid particles, either separately or in mixtures. 
Gaseous pollutants constitute approximately 90% of the total mass 
entering the atmosphere while the other 10% are made up of particulates 
and liquid aerosols [1]. Pollutants emitted directly from identifiable 
sources are classified as primary pollutants. On the other hand, 
pollutants produced by interactions among primary pollutants or by 
reactions with atmospheric constitutents are categorized as secondary 
pollutants [2]. These pollutants can be hazardous to both plant and 
animal life. By and large, the most harmful air pollution is due to 
the production of acid droplets through the interaction of sulfuric and 
nitric oxides with humid air. These oxides evolve primarily from com- 
bustion of fossil fuel and refuse. 

A complete description of polluted air masses is quite difficult 
since generally the polluting entities do not retain their exact 
identities after entering the atmosphere. Pollutants undergo thermal 
and photochemical reactions with continuous spatial and time 
variations. In the vicinity of polluting sources locally high con- 
centrations and significant changes in time may occur. Subsequent dis- 
persion depends strongly upon the prevailing winds, atmospheric stability 
and topographic features. Levels of pollutants are determined in terms 
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of concentration density, i.e., concentration per unit volume. Similar 
health hazards can arise from a short time exposure of a few minutes to 
locally high concentrations or of long time exposures of hours and 
even days to low concentrations. 

Detection of pollution through direct sampling is rather complex 
sines a great number of sampling stations distributed over large areas 
are needed. In the United States there are over 7,000 sampling stations 
[1]. Most of these stations, however, are located in metropolitan and 
industrial areas, i.e., in area sources and in the vicinity of major 
point sources such as high smoke stacks. To monitor the spatial and 
time variations of pollutants to obtain contour maps of trace constituents, 
surveys within relatively large distances downstream of the sources are 
required. For instance, aerial photographs of power plants reveal that 
the most dangerous surface pollution concentrations occur within 5 to 10 
km downwind of the plants in the case of a limited mixing plume [3] . Even 
surface concentrations at distances beyond 30 km from the source were 
detected for fanning and inversion breakup plumes [3], Surveys of plume 
rise and dispersion over such large areas using sampling stations are 
costly and inherently difficult due to the continuous changing meteorological 
conditions. Moreover, the use of sampling stations to monitor line sources 
such as highways and runways (mobile sources), and unpredictable downwind 
precipitation areas is practically precluded due to the extreme range of 
area involved and the wide changes in concentration. In many cases, 
furthermore, the setup of sampling stations is prevented by local safety 
regulations, e.g., erection of tall meteorological towers in vicinity of 
airports . 

Pollutants are deposited into the atmospheric boundary layer where 
they are conveyed either as active or passive scalars by advection and 
turbulent diffusion. In the atmosphere the turbulent diffusion depends 
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on both shear flow and vertical temperature distribution. The buoyancy 
forces, which are determined by the temperature variation, play a prime 
role in the dynamics of air motion. When stable conditions prevail, i.e., 
subadiabatic temperature lapse rates, the turbulent diffusion of con- 
taminants Is retarded by buoyancy forces. The extreme case of stable 
conditions is represented by temperature inversions. Strong and stable 
inversion layers yield extreme circulation conditions where vertical 
turbulent mixing of pollutants is effectively prohibited by the tempera- 
ture lid. Dangerous inversion which can persist for several days 
usually develop within high pressure ranges. Over urban areas tempera- 
ture inversions are commonly found beyond a certain mixing depth [4, 5], 
Generally, the urban inversions arise in the form of a dome-shaped 
envelope due to the uneven heat distribution. Furthermore, on clear 
mornings at sunrise and often or. evenings at twilight with light vinds 
or calm conditions near the surface, a base inversion layer occurs. This 
layer may extend from the ground up to heights of 150 to 750 m depending 
upon the local topographic configuration [6]. Simultaneously, a second 
inversion layer may exist at higher altitudes [6]. Inside an inversion 
layer the highest surface concentrations occur within 5 to 10 km down- 
wind of the sources and can persist for prolonged time depending upon 
the type and lifetime of the inversion [3]. 

Air pollution surveys are strongly dependent on detailed knowledge 
of prevailing turbulent winds and atmospheric stability conditions 
within relevant large areas. Remote-sensing technology possesses the 
potential to monitor both meteorological parameters and needed data on 
the dynamics of mixing of pollutants in sufficient detail and within 
substantial large areas. The orbital stereo photographs taken during 
the Gemini missions clearly substantiate the use of stereo techniques 
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to survey contaminated air masses over exceedingly large areas [7]. 
Furthermore, advance detection of unexpected pollution is feasible. As 
a result, adequate regional air quality standards and control regulations 
can be set forth. However, to determine chemical composition of 
pollutants, size and distribution of particulates and health endangerment 
levels, direct sampling is needed. Remote sensing by means of stereo- 
scopic observations is capable of yielding a host of information which 
can. augment manifold. the knowledge on the dynamics of mixing of pollutants 
and surface concentrations . 

ANALYSIS OF STEREOSCOPIC IMAGES 

Stereoscopic viewing of targets essentially reduces to stereo 
comparison of a pair of images with some overlapping areas . Such two 
successive photographs can be taken by a precision aerial camera along 
the flight path of an aircraft. TWo pictures taken in succession and 
the corresponding overlapping areas indicated by shaded regions I and II 
are shown in Fig. 1(a). Next, two arbitrary objects are selected within 
the area of overlap. They are denoted by 1 and 2 in the first photograph 
while their conjugate images in the second picture are designated by 
1' and 2', respectively. These objects can be a sharp ground feature 
and a tenuous but visible polluted air mass as portrayed in Fig. 1(b). 

Under the assumption of a fixed (stationary) pollution target, the height 
and three-dimensional extent of the plume can be determined by stereoscopic 
comparison of the two images [8]. Titus, the time lapse between two 
successive exposure's must be short enough such that the distance traveled 
by the pollution target can be disregarded. The human matching of the 
stereo images is carried out using a stereocomparator. Subsequent 
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attempts to align the pollution target will lead to a mismatch for the 
ground object. This misalignment is caused by the fact that the two 
targets are situated at different heights. The distance difference, i.e., 
the parallactic distance, can be estimated from the mismatch of the 
ground reference object as illustrated in Fig. IQ}). Then, the height 
of the pollution target is computed by triangulation utilizing the 
parallactic distance and the difference in the view angles from the two 
images for the pollution target (the parallactic angles). This method 
can also be used to estimate the features of an inversion layer. When 
the motion of the polluted mass is not negligible, a third stereo image 
is needed. Then, two independent parallactic distances are obtained. 
Similar triangulations yield both plume altitude and its horizontal 
velocity. 

Results using this stereo analysis to deduce the three-dimensional 
extent of visible plumes are reported in Ref. 9. The three-dimensional 
extent of plumes were estimated using vertical and oblique stereo image 
pairs with 60% overlap taken at altitudes of 10,000 and 15,000 ft, 
respectively. The accuracy of this method decreases drastically with 
distance from the source. At 3 miles downwind of the polluting source, 
the error in approximating the plume height is ± 100 ft [9]. Beyond 
this range, evaluation of plume height becomes uncertain. Human 
photogrammetric stereo comparison is limited even under optimum conditions. 
For barely visible and invisible plumes such as oxide gases, human stereo 
interpretation is completely impractical. With increased use of pre- 
cipitators in high smoke stacks, the visibility problem becomes more 
compound. 
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The stereoscopic analysis of images can be highly improved and refined 
when adequate digital cross -correlation techniques are utilized. Ini- 
tially, these techniques were developed and employed for remote measure- 
ment of sound sources, heat radiation, wind, and turbulence using passive 
optical crossed beams [10, 11, 12]. The crossed-beam correlation method 
was further applied to active sensors. A review of the applications of 
electromagnetic correlation techniques can be found in Ref. 13. 

Basically, this method combines remote measurements with adequate data 
processing and statistical analysis. 

To carry out digital correlation of a stereoscopic image pair, 
conversion of each picture into discrete scan lines is necessary [14]. 

These lines, which are depicted in Fig. 1(a), can be obtained by scanning 
stereo photographs with a microdensitometer [9, 15]. Essentially, the 
discrete scan lines constitute the scan history of the image in the 
form of a digitized record of the radiation power within a certain field 
of view. The scan history consists of a sufficient number of 
adjacent scan lines which supply an image of the target area in terms of 
a strip map. It is, further, important to note that scan histories 
within wavelength ranges where photographic emulsions are not available 
can be supplied by utilizing electromechanical stereoscopic scanners. 

The scan history can be expressed in terms of a signal 

x(t) • x(m,n) , (1) 

where m designates the scan line number relative to the position of 
a ground reference object whose scan line is m * 0. The sample point 
number along a particular scan line is denoted by n . Tt is determined 
with respect to the sample point n * 0 that intersects the ground 
reference feature. 
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The stereo analysis is performed by computing the cross -covariance 1 

function of the two overlapping areas I and II. To calculate the cross - 
covariance, one of the images. *■%., image II, is delayed by k scan 
lines and l sample points. Then, the sample cross -covariance estimate 
is 

. M N 

C (k,t) * 4 M}T l l x(m,n)y(m*k,n+t) , (2) 

** -M -N 

where x(m,n) is the signal from image I, and y(m+k,n+£) is the 
delayed signal from the overlaying area II. The number of resolution 
elements per strip is denoted by 4MN. A direct quantitative measure 
for the matching is provided by the cross -covariance. The best match 
of the features related to the pollution mass is indicated by the peak 
value of the cross -covariance. In a similar manner, but independently, 
the match of the ground reference can be obtained. As a result, the 
parallactic distance is obtained and, then, the plume features (or 
inversion characteristics) are evaluated by triangulation as portrayed in 
Fig. 1(b). 

The signals from stereo images are random variables of time and, 
thus, they can be expressed as ! 

x(t) « Xj (t) ♦ x 2 (t) , (3.1) j 

and 

y(t) - >’j(t) ♦ y 2 (t) , (3.2) 

where the subscripts 1 and 2 denote the target signal and transmission 
noise, respectively. Then, the cross -covariance is 

C xy (T) ■ XjtO y^t ♦ t) ♦ Xjft) y 2 (t ♦ t) 

♦ x 2 (t) y x (t ♦ t) ♦ x 2 (t) y 2 (t * tT . (4) 


i 
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where the t designates the time delay, and the overoar denotes time 
integration over the entire scan history. The first term on the right- 
hand side of Eq. (4) is the cross -covariance of the common signals. It 
is important to remark that this product is always positive since botn 
signals are of same sign being generated by the same random process 
[12]. In other words, both stereo images are affected in a similar way 
by the polluted air mass target. This is true provided that t!'« motion 
of the pollution target is disregarded, i.e., a stationary target. 

The other three products represent cross -covr. nances of uncorrelated 
random variables caused by statistically independent random processes. 
Consequently, they vanish provided that the integration time is Ic-.g 
enough. Basically, the cross -covariance is dominated by the product 
^ . This term increases linearly with integration time. On the 
other hand, the other three terms produce oscillations with equal 
likelihood about the average value, i.e., negative and positive oscilla- 
tions. Their contributions to the time integral increase with the 
square root of integration time [12]. Adequate averaging methods to 
cancel the noise by product integration for finite integration time are 
proposed in Ref. 12. 

The mean square error of the cross -covariance estimate due to 
indiv dual resolution element (or finite integration time) is given by 
the variance, 

v * r l c xy (M)] ■ 4 j^rj- l [x(m,n)y(m*k,n*t) - C xy ] 2 , (5) 

S|A 

since C (k,4)] is an unbiased estimate of the true value ot the 

cross -covariance (T [16]. Moreover, the mean square error decreases 

*/ 

with increasing integration time. Next, the confidence to obtain a 
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meaningful match is expressed by the ratio of the cross-covariance 
estimate to its root-mean-square (rms) error 


a 


k,£ 


✓Var [C xy Ck,£)} ’ 


(6) 


which is the coefficient of confidence. 

In practice, additional limitations exist due to the uncontrolled 
effects of the optical environment. Averaging techniques employing the 
variations of piecewise mean values to recognize and eliminate these 
effects were attempted [12]. A sample of the instantaneous product xy of 
two total signals for two different wavelength ranges from a stereo scanner 
is shown in Fig. 2. The continuous strongly oscillating curve represents 
the scan time history. The cumulative average of the instantaneous product 

i T 

C (k,£,T) = ± f x(M,t) y(k,A,t)dt , (7) 

*y 0 

is a function of the variable record length (or integration time) , 

0 < T < 2MN. With increasing integration time the variations of the 
cumulative average diminish since the changes in the instantaneous 
product cancel each other. The negative fluctuations of the instantaneous 
product represent signal components which are unrelated to the target 
since the contribution of the common signals from the polluted mass to 
the instantaneous product is always positive. Moreover, the positive 
instantaneous products also contain signals which are uncorrelated to 
the pollution target. These positive products cancel the negative 
products. The shaded areas in Fig. 2 delineates the power of the 
unrelevant signals. Their cancellation through the integration process 
is indicated by the decreasing change in the variations of the cumulative 
average. Furthermore, with increasing number of instantaneous products 




i 

l 
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the mutual cancellation of the unrelated signals is enhanced. This 
process of cancellation by product integration is comp’eted when the 
cumulative average reaches a constant value [12, 14] Practically, this 
value represents the power of the common signals from the pollution 
target. Consequently, the product integration process permits recovery 
of very weak signals provided that the target is sufficiently large 
and homogeneous. The latter condition is needed to allow suppression 
of the noise [12]. This method can eventually be further utilized for 
a nonhomogeneous polluted mass by subdividing it into homogeneous domains 
using unsupervised classification techniques [14], As a result, even 
tenuous pollution targets are eventually retrievable using this correla- 
tion analysis. Furthermore, this technique can be applied to evaluate 
the water vapor burden within a column of air. 

EXPERIMENTAL RESULTS 

To demonstrate the feasibility of stereoscopic remote sensing of 
polluted air masses and of stereo-correlation analysis of stereoscopi 
images, an extensive field experimental program ir- currently being 
conducted. To start with, it was necessary to ascertain the capability 
of using stereo images to detect the planar extent of a moving pollution 
target. Subsequently, it is planned, during the second stage of the 
program, to assess the three-dimensional features of a moving pollution 
target utilizing the correlation analysis. The objectives of a 
systematic field experiment could not be achieved by simply surveying 
pollutants randomly emitted into the atmospnere. It was imperative to 
surmount the unpredictable effects of continuous meteorological condi- 
tions. Furthermore, it was desired to inject pollutants under con- 
trolled circumstances into a known flow environment and unchanged flow 
conditions over a long enough time period. 
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For these reasons, it was decided to simulate the environmental 
conditions using the wake flow produced by a 10 ft diameter fan of 
variable pitch (Hartzell A120-6). This fan constitutes the core of 
the Colorado State University Environmental Field Station. The fan, 
which is driven by an internal combustion engine, could produce 
continuous varying air speeds up to 20 mph. Within the wake large scale 
tumulence is produced and, consequently, the environmental conditions 
are adequately simulated for pollutant injection. To reduce undesirable 
interferences with the ambient winds, all experiments were conducted 
under calm conditions or light winds up to 2 mph at the most. Continuous 
monitoring of the prevailing winds was performed using cup anemometers. 

Contaminants in the form of gases, liquid and solid particulates 
can be released under governed conditions from point sources and/or line 
sources. The pollutants may be emitted either separately or premixed in 
any desired ratios. Locations of polluting sources were selected based 
on velocity distribution and turbulence characteristics within the 
wake to avoid any perturbations generated by the fan supporting 
structure. The sources, consequently, were located at least 20 ft 
downstream of the fan. Moreover, since the sources are mobile they can 
be located in any configuration and at any position of interest. It 
is worth pointing out that similar sources can be utilized to generate 
controlled water vapor plumes within the wake. 

The Environmental Field Station is located on a flat surface free 
of immediate natural or artificial obstructions. At the field station 
almost daily base temperature inversions occur on clear mornings and 
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in the evenings at twilight. The station is equipped with an Analog 
Data Acquisition System. In addition, a Digital Data Acquisition 
System is available at the site for on-line data reduction and analysis. 
Fixed ground reference targets needed for the stereo photographs were 
set up. 

To ascertain the practical applications of the stereoscopic remote 
sensing, corresponding ground measurements are necessary. Furthermore, 
comparison of the results will lead to development of adequate calibration 
techniques for the stereo remote sensing. Ground measurements of 
velocity and turbulence were carried out utilizing a novel three- lead 
hot-wire anemometer system conceived, designed and built at Colorado 
State University [17]. Generally, the use and performances of hot-wire 
anemometers depend on the length of the cables connecting the hot wire 
(the sensor) to the anemometer's bridge. The length of the connecting 
leads and, hence, their resistance affects strongly the bridge balance, 
the noise level and the frequency response of the anemometer. The accepted 
length of the leads is limited to about 25 ft. Whenever longer cables 
are utilized adequate compensation is required to allow the bridge 
balance. However, both the noise level and frequency response are 
increasing and decreasing, respectively, with the cable's length. The 
new three-lead system removes the cable's length constraint since 
connecting leads in excess of 500 ft length can be utilized without 
need for compensation [17]. In this system the balance of the 
bridge is not dependent upon the length of the connecting cables. The 
noise level remains practically unaffected by the lead's length, and 
the decrease in frequency response is negligible. The new remote system 
affords innumerable applications of hot-wire anemometry techniques in 



atmospheric measurements. Specific hot-wire interpretation methods 
developed for flow with large fluctuations, such as encountered in the 
atmosphere, are employed [18]. 

The aerial stereoscopic photographs were taken employing a Wild 
RC8 automatic precision mapping camera system with 93° coverage (Wild- 
Heerbrugg Inst.). The camera is equipped with a 6 in. Universal Aviogon 
lens cone (f/S.6) with rotary shutter. An antivignetting filter Wild 
A.V. 2.2 x (Wild-Heerbrugg Inst.) is utilized to provide uniform lighting 
across the lens and, thus, uniform exposure on the film is obtained. 

All pictures were taken using Kodak aerocolor negative estar base film. 
Included in the camera system are a driving mechanism, a control box, 
a viewfinder telescope, a clock, and an altimeter. A general view of the 
camera system is shown in Fig. 3. The camera produces 9 x 9 in. photo- 
graphs with remarkably low distortion and high resolution. The former is 
smaller than 0.01 mm whereas the latter varies from 50 lines/mm in 
the center to 25 lines/mm in the comers. This high resolution 
permits adequate conversion of the picture into discrete scan lines for 
further digital correlation analysis. The film advance per exposure 
amounts to 10 in. allowing readings of the incorporated clock and alti- 
meter to be recorded on each photograph. Continuous adjustment of the 
exposure time from 1/1000 to 1/7000 sec is provided by the driving 
motor. Four overlap areas of 20, 60, 70 and 80%, respectively, are 
available. The shortest time interval between exposures, i.e., the 
time lapse, is roughly 3 to 3.5 sec. A plume conveyed by such a 
velocity that the distance traveled during the time lapse is small 
compared with its overall horizontal extent will appear on two 
successive photographs as being approximately stationary. At this 







14 


point, it is important to notice that the plume is "stationary" with 
respect to stereoscopic photography since it will roughly affect a 
stereo-image pair in the same manner as a ground reference object. Thus, 
one can consider the plume as being stereoscopical ly stationary. Such 
situations occur particularly within inversion layers, light winds or 
calm conditions. For aerial operation, the camera system is flown in 
Colorado State University Aero-Commander 500B research aircraft. 

Direct measurement and recording of the space coordinates of 
any desired image on a stereo-photograph pair is carried out utilizing 
the Wild STK-1 Stereocomparator/Digitizer/IBM Card Punch (Wild-Heerbrugg 
Inst.). The coordinates of selected objects can be measured with a 
precision of lym (1/25,400 in.). It is important to remark that the 
evaluation of the vertical coordinates may be arduous due to the 
haze of the plume. Scanning by means of a microdensitometer, which 
measures the transmittance of the photographic emulsion, can undoubtedly 
supply more reliable estimation of the vertical extent of the pollution 
target. Such a scanning is being planned. 

Henceforth, the preliminary results of the first stage of the 
feasibility investigation are presented. To begin with, intensive, 
although not exhaustive, visualization studies of circulation of 
visible plumes generated by colored smoke were conducted. The smoke 
point sources were located at various positions along the streamwise 
direction of the wake flow and at several elevations above the ground. 
Several frames from a movie showing the smoke plume circulation, en- 
trainment and dispersion without and with temperature inversion are 
given by Fig. 4. The fan is also shown in Fig. 4(a). Relatively 
large scale vortices were clearly discerned when unstable conditions 
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(superadiabatic lapse rates) prevailed. The photographs shown in Figs. 
4(a) and 4(b) were taken under such conditions. Without temperature 
inversion, the smoke plume is dispersed quickly at higher altitudes due 
to the vortices within the wake. However, in the neighborhood of the 
sources, due to the large scale vortices, deposit of pollutants on the 
ground can occur as indicated by the smoke circulation observed in 
Fig. 4(b). Under inversion conditions small scale vortices were 
prominent as illustrated by the photographs shown in Fig. 4(c). Further- 
more, relatively high smoke concentrations persisted for a prolonged time. 
Knowledge of the scale of vortices is essential for the study of both 
dispersion and concentration of pollutants, and temperature distribu- 
tion. Within the earth -atmosphere interface, in particular, the 
prevailing vortices play a prime role in the surface concentration of 
pollutants. 

Next, numerous overflight surveys of the smoke plume were carried 
out. Stereoscopic photographs were taken from several altitudes. All 
photographs were obtained with 60% overlap and with a time lapse of 3 
sec. The aerial photogrammetric surveys were conducted either under calm 
conditions or light winds up to about 1 mph and temperature inversion. 

A stereo-image triplet of smoke circulation taken from an altitude of 
2000 ft is shown in Fig. 5. Analysis of the stereo images by means of a 
stereocomparator revealed that the maximum measurable windwise extent of 
the plume was about 70 ft. The shape and extent of the smoke plume in the 
horizontal plane is depicted in Fig. 6. For a wind speed of 1.5 ft/sec 
(« 1 mph) , the smoke plume traveled a distance of 4.5 ft during the 
time lapse. Since this distance is about 6.4% of its windwise range, 
one can conceivably assume that, for most practical purposes, the plume 
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could be considered stereoscopical ly stationary. The similar extent of 
the smoke plume obtained by comparison of a series of stereo images 
taken under same conditions substantiates the foregoing conclusion. 

Basically, the smoke plume was produced by injecting smoke into the 
wake. The smoke was then transported by the flow and practically entrained 
throughout almost the entire wake. Hence, the wake can be considered as 
being an invisible plume whose overall spatial extent roughly 
coincides with that of the visible smoke plume. To corroborate the 
horizontal extent of the smoke plume deduced from the stereo images, 
hot-wire anemometer ground measurements of the planar extent of the 
corresponding invisible plume were performed. These measurements 
were conducted under same wind and inversion conditions as for the 
remote stereo survey but not simultaneously. The boundary of the 
invisible plume, i.e., the wake boundary, was defined as the location 
where the local streamwise velocity reduces to about 2% of the fan rotor 
tip velocity [19]. The results of the hot-w'ire anemometer survey are 
shown in Fig. 6 togethei with the horizontal extent of the visible smoke 
plume. A striking similarity between the extents of the visible and in- 
visible plumes is clearly discerned. Their planar extents coincide 
within ± 10%. This consistent general congruence substantiates the 
feasibility of using stereo images. Note that this agreement was 
obtained for a relatively small, tenuous plume as clearly indicated 
by Fig. 5. Thus, stereo images can be used to survey even barely 
visible and weak plumes. Currently, the evaluation of the three-dimen- 
sional features of the plume, i.e., stereoscopic estimation of the height 
using digital correlation analysis, is underway. 
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CONCLUDING REMARKS 

The feasibility of stereoscopic remote sensing of polluted air masses 
by means of flown stereo cameras is indicated by the preliminary results 
presented in this work. The evaluation of the planar extent of a plume 
based on both stereo images and ground measurements exhibits a reason- 
able agreement. Even tenuous plumes of relatively limited spatial extent 
can be monitored using stereo images. In this experiment the windwise 
extent of the plume was roughly 70 ft. 

Analysis of a stereo-image pair can be used to establish the features 
of a moving polluted air mass provided that the condition of stereoscopic 
stationarity is satisfied. When the distance traveled by a plume 
during the time lapse between two successive exposures is small compared 
with its overall horizontal extent, one can assume *hat the plume is 
stereoscopical ly stationary. This situation usually occurs within tempera- 
ture inversion layer, under light winds or calm conditions. When the 
condition of stereoscopic stationarity is not fulfilled, e.g., strong winds, 
a stereo-image triplet may be necessary. Currently, analysis of such 
triplets is being attempted. 

By and large, standard stereo comparison analysis is limited due to 
variations in the concentration of pollutants within the plumes, the 
haze of the plumes, and the continuous changing optical background. 

Digital stereoscopic analysis based on correlation techniques may possess 
the potential to overcome these difficulties. An intensive effort to 
demonstrate the efficacy of digital stereo correlation is planned. It is 
expected that the completion of the correlation analysis of stereo-image 
pair and triplet, which is currently in progress, will hopefully provide 
the needed support for the feasibility of remote stereoscopic survey of 
polluted air masses. The development of mathematical models for long 


term pollution distribution ar.C '•.onrentraiit i over large areas depends, 
to a large extent, un r’.ttai i*d knnwis.lgt of both meteorological condi- 
tions and polluted air mass; feat<x:cs' . Stereoscopic techniques eventually 
can be used efficiently tc pr-v iJe contour maps of pollution including 
isoconcentration density curves in addition to meteorological parameters. 
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Fig. 1 (a) Stereoscopic imagery, (b) Determination of 

pollution target height by stereoscopic comparison 
and correlation techniques. 
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Fig. 2. Cancellation of unrelated signals by the 
product integration method. 
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Fig. 4 Visualization of smoke plume circulation under unstable conditions: (a) point 

source upstream of the fan; and, (b) point source at 40 ft downstream of the fan; 
and within a temperature inversion layer: (c) point source at 150 ft downstream 
of the fan. 






Fig . 5 Reproduction of colored stereo- image triplets of smoke 
circulation talon from an altitude of 2000 ft. 




